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1.  INTRODUCTION 


One  of  the  goals  of  the  We^ns  Technology  Directorate  of  the  U.S.  Army  Research  Laboratory  is 
to  understand  the  penetration  performance  of  penetrators  under  certain  circumstances.  The  penetration 
depth  of  tungsten  rods  into  rolled  homogeneous  armor  (RHA)  is  one  such  performance  measure.  It  was 
shown  by  Frank  and  Zook  (1987)  that  there  exists  an  impact  velocity  such  that,  for  a  constant  energy, 
constant  length-to-diameter  (L/D)  ratio  rod,  the  penetration  depth  is  a  maximum.  This  optimum  velocity 
can  be  obtained  from  the  penetration-vs.-velocity  curve  in  a  simple  manner  (Frank  and  Zook  1987, 1991). 
This  has  been  done  for  three  different  analytical  representations  of  data  for  tungsten  penetrators  (de  Rosset 
1992).  All  three  data  representations  result  in  an  optimum  velocity  of  about  2.1  km/s.  However,  none 
of  the  three  representations  is  well  documented  with  respect  to  the  data  used  or  the  method  by  which  the 
constants  for  the  analytical  expressions  were  derived.  In  addition,  only  one  of  the  analytical  expressions 
contains  dependence  of  penetration  on  length  and  diameter.  The  puipose  of  this  woik  is  to  examine 
several  functions  which  might  be  used  to  fit  penetration-velocity  data  for  tungsten  penetrators  impacting 
RHA.  Each  of  the  functions  will  have  an  explicit  dependence  on  L/D  ratio.  The  database  will  be 
screened  for  acceptable  data  (primarily  low-yaw  tests)  and  will  feature  a  range  of  L/D  ratios.  These 
functions  will  then  be  used  to  determine  an  optimum  velocity,  based  on  the  fitting  parameters  for  each 
function.  The  dependence  of  the  optimum  velocity  on  L/D  ratio  will  then  be  apparent. 

The  approach  used  was  first  to  gather  and  screen  data.  One  of  the  primary  data  sources  used  for  this 
woiic  was  a  recently  published  compendium  of  terminal  ballistics  data  (Anderson,  Morris,  and  Littlefield 
1992).  Other  reports  containing  tungsten-vs.-RHA  data  were  also  examined  (Magness  and  Leonard  1993; 
Frank  and  Zook  1987, 1990, 1991, 1992).  Next,  a  suitable  statistical  analysis  program  was  selected.  The 
work  started  out  using  RS/1,  which  was  soon  discarded  in  favor  of  the  more  user-friendly  SigmaPlot. 
Finally,  several  fimctions  were  considered  to  represent  the  data.  The  results  from  the  fits  to  the  data  are 
compared  in  the  discussion  sectioa 

Section  2  describes  the  general  features  of  SigmaPlot,  gives  the  criterion  by  which  data  were  accepted, 
and  documents  the  data  used  for  this  study.  The  following  section  describes  the  different  fitting  functions 
considered  in  this  study.  After  the  results  section,  the  findings  of  this  study  are  discussed.  A  concluding 
section  summarizes  the  report 
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2.  APPROACH 


SigmaPlot  (Norby  et  al.  1986)  is  a  program  that  uses  a  wide  variety  of  graphing  techniques  to  fit  data 
in  a  statistical  manner.  The  data  worksheet  in  SigmaPlot  consists  of  over  64,000  rows  and  16,000 
columns.  The  following  parameters  are  used  as  input  to  fit  the  prescribed  curve:  the  ratio  of  penetration 
depth  to  initial  rod  length  (P/L),  penetrator  material  density,  rod  L/D  ratio,  and  rod  striking  velocity  (v). 
The  target  density  was  taken  to  be  7.85  g/cm^  in  all  cases.  The  prescribed  function  is  then  entered  into 
the  program.  In  aU  cases,  the  function  is  limited  to  two  fitting  parameters,  which  are  labeled  A  and  B. 
The  Marquardt-Levenberg  algorithm  is  then  used  by  the  program  to  find  the  values  of  A  and  B,  which 
minimize  the  difference  between  the  data  and  the  function  in  a  least-squares  manner.  The  program  also 
calculates  the  standard  error  for  A  and  B.  This  is  used  to  detemiine  the  coefficient  of  variation  (CV), 
defined  by 

CV  =  standard  error  x  100/parameter  value.  (1) 

CV  is  a  measure  of  how  weU  the  choice  of  the  function  and  the  fitting  parameters  represent  the  data  and 
wiU  be  used  to  select  the  best  fitting  function  of  the  ones  examined. 

All  the  data  were  screened  to  eliminate  high-yaw  impacts.  This  was  done  in  order  to  reduce  the 
amount  of  scatter  and  assure  that  the  maximum  penetration  at  each  velocity  was  being  entered  into  the 
database.  The  criterion  used  to  exclude  data  was  that  the  impact  yaw  could  not  exceed  a  critical  yaw 
given  by  the  following  formula  (Bjerke  et  al.  1992): 

critical  yaw  =  sin"^  [(H-D)/2l].  (2) 

Here,  H  is  the  penetration  channel  entrance  diameter,  L  is  the  penetrator  length,  and  D  is  the  penetrator 
diameter.  The  channel  entrance  diameter  is  also  found  in  Bjeike  et  al.  (1992)  and  given  by 

H  =  (l.524  +  .3388v  +  .1286v2)d,  (3) 


where  v  is  the  impact  velocity. 
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An  attempt  was  made  to  screen  the  targets  according  to  their  hardness.  This  was  not  always  possible, 
since  in  some  instances  the  target  hardness  was  not  reported.  In  those  tests  where  the  hardness  was 
reported,  the  values  ranged  between  BHN  255  and  290.  Targets  identified  as  mild  steel  were  not  included 
in  the  database. 

The  data  points  used  in  this  study  are  provided  in  the  Appendix  in  the  event  that  others  may  wish  to 
use  them  to  fit  different  functions.  Table  1  summarizes  the  sources,  range  of  parameters,  and  amount  of 
data  from  each  of  the  sources.  The  data  set  is  limited  to  0°  obliquity  targets,  and  the  range  of  penetrator 
masses  was  finom  30  to  250  g.  While  this  data  set  does  not  include  every  test  tiiat  has  been  conducted, 
it  does  represent  a  large  number  of  tests. 


Table  1.  Data  Summary 


Source 

L/D 

Velocity  Range 
(m/s) 

No.  of  Data  Points 

Hohler  and  Ship  (1991a) 

1,  10,  16.3,  17.5,  20, 
21.7,  22,  22.5,  32 

951-3,663 

55 

Holder  and  Stilp  (1991b) 

1,  3,  5,  7 

2,281-3,652 

18 

de  Rosset  et  al.  (1989) 

15,  20,  30 

2,140-3,050 

13 

Magness  and  Leonard 
(1993) 

20 

1,167-1,680 

8 

Sfisby  (1984) 

23 

1,865-^,525 

10 

Frank  and  Zook  (1990) 

1.0 

906-4,881 

18 

Figure  1  shows  this  data  graphically. 

3.  DATA  REPRESENTATION 

In  choosing  a  function  to  fit  the  penetration  depth  vs.  velocity  for  tungsten  impacting  RHA  at  0° 
obliquity,  there  were  several  criteria  which  had  to  be  met  First,  it  was  desirable  that  the  function  had 
some  connection  to  physical  reality.  That  is,  a  high-order  polynomial  was  not  considered.  (See  Bless 
et  al.  [1994]  for  an  eight-parameter  fit.)  However,  for  this  particular  study  the  effect  of  target  and 
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Figure  1.  Graphical  representation  of  penetration  data. 


penetrator  hardnesses  was  ignored.  This  effect  would  have  to  be  included  in  future  refinements  of  the 
model,  assuming  that  there  is  enough  data  in  existence  where  the  target  hardness  was  measured  (and  the 
target,  in  fact,  had  a  uniform  through-thickness  hardness).  The  function  had  to  contain  a  dependence  on 
L/D.  The  behavior  of  the  function  at  input  variable  extremes  (e.g.,  very  high  velocity,  very  high  or  low 
L/D  ratio,  etc.)  had  to  produce  the  expected  limiting  conditions.  The  number  of  fitting  parameters  was 
limited  to  two.  Finally,  the  function  had  to  fit  the  data  in  a  statistical  sense. 


The  first  function  which  was  examined  was  suggested  by  Frank  (1994): 


P/L  =  l/p»  (v/A)f^'^^^  exp  [-(B/v)^ 


(4) 


where  p  is  the  square  root  of  the  ratio  of  target  and  penetrator  densities.  In  this  and  the  functions  which 
follow,  p  is  present  to  account  for  small  differences  in  various  tungsten  alloy  densities.  It  also  shows 
explicitly  that  the  formulation  approaches  the  density  law  for  high  L/D  ratios  and  high  velocity.  This 
function  produces  the  characteristic  v^^  dependence  for  L/D  =  1  penetrators  at  high  velocity,  and,  the 
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expected  dependence  for  L/D  0  (Orphal  et  al.  1993).  The  exponential  is  that  used  by  Lanz  and 
Odermatt  (1992)  to  fit  full-scale  data.  A  more  thorough  discussion  of  the  exponential  function  is  given 
by  Frank  et  al.  (1992),  who  relate  B  to  such  parameters  as  target  resistance  and  penetrator  strength. 
Equation  4  wiU  be  referred  to  as  Formula  1. 

The  special  feature  of  this  fimction  is  that  a  simple  analytical  expression  for  the  optimum  velocity  of 
constant  energy  penetrators  can  be  obtained.  Assuming  constant  penetrator  energy,  constant  penetrator 
geometry  (L/D  ratio),  and  penetration  depth  equal  to  penetrator  length  times  a  function  of  penetrator 
velocity  and  other  relevent  parameters  (density,  material  strength,  etc.),  Frank  and  Zook  (1991)  give  the 
condition  for  optimum  penetrator  velocity  as 

v(aP/av)  =  2/3  P.  (5) 

For  Formula  1,  the  optimum  velocity  is  given  by 


'opt 


=  B  /T" /v[i-(3/(1+2L/D)))  , 


(6) 


which,  for  L/D  — >  °o,  is  identical  to  the  optimum  velocity  derived  by  Frank  et  al.  (1992). 

Formula  1  was  used  to  fit  the  data,  and  it  was  found  that  the  fit  at  low  to  moderate  velocities  was  fair. 
A  second  formulation  was  examined  in  an  attempt  to  improve  the  fit.  This  formulation  divided  the 
penetration  into  two  parts:  a  steady-state,  long-rod  contribution  and  a  low  L/D  portion.  The  formulation 
follows  the  concept  of  Christman  and  Gehring  (1966),  who  divided  the  penetration  process  into  four  parts. 
The  two  parts  of  the  current  formulation  correspond  to  the  primary  and  secondary  penetration  phases  of 
Christman  and  Gehring.  The  formulation  is  given  as 


P/L  =  (1/p) 


((1  -D/L)  exp 


+  D/L 


(7) 


The  idea  here  is  that  a  length  of  rod  equal  to  L  -  D  contributes  to  and  is  eroded  in  the  steady-state  phase, 
while  the  remaining  length  D  contributes  to  the  final  portion.  For  L  =  D,  the  usual  dependence  of  P  on 
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is  recovered.  The  approach  is  similar  to  that  used  by  Charters  and  Oiphal  (1988)  and  is  denoted  by 
Formula  2. 

This  representation  results  in  an  expression  for  the  optimum  velocity  which  is  even  simpler  than  that 
given  by  equation  6.  In  this  particular  case, 

Vopt  =  (8) 

Thus,  while  P/L  depends  explicitly  on  L/D,  the  optimum  velocity  does  not. 

Charters  (1992)  has  concluded  that  a  function  of  the  form  (v  -  vo)^^  fits  L/D  =  1  penetration  data 
better  than  a  simple  v^^  form.  Frank  and  Zook  (1990)  have  also  derived  an  analytic  expression  for  low 
L/D  penetrators  which  fits  the  data  well  at  aU  velocities.  The  general  velocity  dependence  of  their  model 
was  used  to  replace  the  v^^  in  the  previous  model,  now  denoted  as  Formula  3: 

P/L  =  1/p  ((1  -D/L)  exp  [-(B/V)^]  +  D/L  *  In  [l  +(v/A)^]).  (9) 

The  expression  for  the  optimum  velocity  using  this  foimulation  is  not  straightforward.  While  it  may  be 
possible  to  solve  for  v^pj  in  an  analytical  expression,  the  approach  used  here  was  simply  to  obtain  v  opt 
numerically  for  the  desired  L/D  ratio. 

4.  RESULTS 

The  calculated  values  of  A  and  B  are  shown  in  Table  2  along  with  the  standard  error  and  CV. 


Table  2.  Curve-Fitting  Parameters 


The  formulas  and  associated  values  of  A  and  B  were  used  to  generate  the  optimum  velocities.  Table  3 
compares  the  optimum  velocities  for  each  of  the  formulas  at  selected  L/D  ratios. 


Table  3.  Optimum  Velocities 


|R99 

mSm 

■SH 

L/D  10 
(m/s) 

L/D  30 
(m/s) 

Formula  1 

2,516 

2,230 

2,054 

1,950 

Formula  2 

2,104 

2,104 

2,104 

2,104 

Formula  3 

2,702 

2,397 

2,220 

2,120 

The  dependence  of  the  optimum  velocity  on  L/D  for  each  of  the  formulas  is  shown  in  Figure  2.  For 
Formula  1  and  Formula  3,  there  is  a  substantial  change  in  the  optimum  velocities  as  the  L/D  ratio  goes 
from  1  to  5.  For  values  of  L/D  over  15,  there  is  not  a  substantial  change  in  the  optimum  velocity.  The 
optimum  velocity  for  Formula  2  is  independent  of  the  L/D,  as  previously  stated.  The  value  of  the 
optimum  velocity  for  Formula  2  is  close  to  that  obtained  for  high  L/D  ratios  using  Formula  3. 

5.  DISCUSSION 

Each  of  the  three  functions  chosen  to  represent  the  data  provides  a  reasonable  fit,  considering  that  only 
two  parameters  were  used.  In  addition,  a  wide  range  of  velocities,  L/D  ratios,  and  data  sources  was  used 
as  input  This  could  also  account  for  values  of  CV  being  higher  than  desired.  The  three  functions  give 
an  explicit  dependence  of  P/L  on  the  L/D  ratio,  and  this  translates  into  a  dependence  of  the  optimum 
velocity  on  L/D  ratio  for  two  of  the  functions.  From  the  results  generated,  it  is  clear  that  the  optimum 
velocity  does  not  change  a  great  deal  for  rods  with  L/D  greater  than  15.  In  addition,  the  results  indicate 
that  the  optimum  velocity  for  long  rods  found  during  this  study  is  not  a  great  deal  different  from  the 
values  previously  determined  (de  Rosset  1992;  Frank  and  Zook  1987,  1991). 

Conversely,  the  optimum  velocity  for  short  rods  may  be  considerably  higher  than  2.1  km/s.  This 
implies  that  segmented  rods  (with  low  aspect  ratio  segments)  must  be  fired  at  much  higher  velocities  than 
conventional  long  rods  to  achieve  their  optimum  performance. 
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17D 

Figure  2.  Optimum  velocity  vs.  L/D  ratio  for  three  formulations. 

Several  different  subsets  of  the  data  presented  in  the  Appendix  were  examined  with  SigmaPlot  to  see 
if  inclusion  of  certain  data  was  producing  abnoimal  results.  The  L/D  =  1  data  points  were  removed,  with 
little  effect  on  the  results.  In  fact,  the  value  of  CV  for  the  functions  increased,  since  the  standard  error 
for  these  points  was  lower  than  average.  In  another  instance,  only  the  L/D  =  20  points  were  considered. 
While  the  value  of  A  and  B  for  Formula  3  changed  somewhat,  the  value  of  CV  did  not  decrease.  A  plot 
of  the  curve  fit  to  these  data  points  is  shown  in  Figure  3.  Data  points  at  1,500  m/s  and  below  faU  below 
the  fitted  curve.  This  may  be  due  to  the  incorrect  assumption  that  an  amoimt  equal  to  L-D  of  the  rod 
length  is  entirely  eroded  during  the  steady-state  process  at  all  velocities.  In  fact,  at  low  velocities,  only 
a  portion  of  the  rod  is  eroded.  Thus,  the  contribution  to  the  penetration  depth  at  low  velocity  would  be 
lower  than  the  formula  predicts.  If  the  formula  were  fit  only  to  the  high  velocity  data,  the  discrepancy 
would  be  even  more  obvious. 
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1.5 


Velocity  (m/s) 

Figure  3.  Comparison  of  L/D  20  data  with  Formula  3  representation. 

In  Formula  2,  the  dependence  of  penetration  depth  on  for  L/D  =  1  penetrators  also  has  the  same 
problem.  The  relationship  can  be  derived  by  setting  the  penetrator’s  kinetic  energy  (KE)  proportional  to 
the  crater  volume  and  assuming  a  hemispherical  crater  shape.  This  may  have  some  vahdity  at  high 
velocity,  but  at  low  velocity  the  relationship  breaks  down.  Note  that  using  this  proportionality  leads  to 
a  constant  penetration  depth  for  a  constant  energy  L/D  =  1  penetrator  at  any  velocity.  Thus,  using  the  v^^ 
dependency  for  the  L/D  =  1  rods  in  Formula  2  will  not  influence  the  determination  of  the  optimum 
velocity. 

An  obvious  point  to  make  about  this  curve-fitting  exercise  is  that  it  is  valid  only  for  the  range  of  data 
employed.  Thus,  for  L/D  less  than  1,  there  may  be  important  features  that  are  not  captured  by  the 
formulas  used  here.  In  fact,  Bjeike  et.  al  (1991)  have  shown  that  P/L  increases  with  decreasing  L/D,  but 
only  to  the  point  where  L/D  is  approximately  1/8.  For  smaller  values  of  L/D,  P/L  acmally  decreases. 
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Consequently,  the  extemely  high  values  of  the  optimum  velocity  derived  from  Formulas  1  and  3  for 
penetrators  with  L/D  near  1  are  suspect. 

6.  SUMMARY 


The  effect  of  the  L/D  ratio  on  the  optimum  penetration  velocity  has  been  examined  for  tungsten 
penetrators  impacting  rolled  homogeneous  armor  at  0°  obliquity.  Three  different  formulas  were  used  to 
represent  about  125  data  points  obtained  from  the  open  literature,  spanning  a  large  range  of  velocities  and 
L/D  ratios.  The  formulas  were  chosen  based  on  some  radimentary  physical  considerations  and  simplicity. 
Only  two  adjustable  parameters  were  used. 

For  two  of  the  formulas,  the  optimum  velocity  decreased  rapidly  as  the  L/D  ratio  went  from  1  to  5. 
The  other  formula  did  not  show  any  dependence  of  the  optimum  velocity  on  L/D  ratio.  AU  formulas 
indicated  that  for  L/D  greater  than  15,  there  was  very  little  change  in  the  optimum  velocity,  which  was 
calculated  to  be  between  1,900  and  2,200  m/s.  This  range  of  values  is  consistent  with  other  values  that 
have  been  calculated  previously  for  long-rod  tungsten  penetrators. 

If  the  optimum  velocity  dependence  on  L/D  is  derived  from  either  Formula  1  or  Formula  3,  then  the 
optimum  velocity  for  a  segmented  rod  with  L/D  1  segments  will  be  in  excess  of  3,000  m/s  for  the 
conditions  examined,  assuming  no  segment  interactions. 
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APPENDIX:  BALLISTIC  DATA 
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Table  A-1.  Data  from  Holder  and  Stilp  (1991a) 
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Table  A-1.  Data  From  Hohler  and  Stilp  (1991a)  (continued) 
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Table  A-2.  Data  From  Hohler  and  Stilp  (1991b) 


Test  No. 

P  2 

(G/cm^) 

L 

(mm) 

L/D 

P 

(mm) 

P/L 

V 

(m/s) 

6536 

17.6 

5.8 

1 

B 

1.897 

2,281 

6537 

17.6 

5.8 

1 

b9 

1.931 

2,291 

6545 

17.6 

5.8 

1 

mm 

2.276 

2,602 

6524 

17.6 

5.8 

1 

mm 

2.483 

3,459 

6517 

17.6 

5.8 

1 

14.8 

2.552 

3,480 

6516 

17.6 

5.8 

1 

15.0 

2.586 

3,652 

6539 

Hi 

3 

29.5 

1.695 

2,400 

6559 

msSM 

3 

28.6 

1.644 

2,513 

6576 

17.6 

3 

31.2 

mSm 

3,448 

6569 

17.6 

17.4 

3 

32.6 

lii 

3,555 

6543 

■a 

29.0 

5 

B 

2,494 

6540 

■H 

29.0 

5 

1.579 

2,803 

6573 

17.6 

29.0 

5 

mm 

1.648 

3,128 

6575 

17.6 

29.0 

5 

49.8 

1.717 

3,513 

6577 

WM 

IB 

B 

B 

B 

2,553 

6556 

■n 

11^ 

2,511 

6508 

■9 

29.4 

7 

49.4 

1.680 

3,257 

6621 

■B 

29.4 

7 

48.6 

1.653 

3,494 
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Table  A-3.  Data  From  de  Rosset  et  al.  (1989) 
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Table  A-5.  Data  From  Silsby  (1984) 


Test  No. 


(G/cm^) 


P 

(mm) 


155.83 

121.79 

121.79 

155.83 

155.83 

121.79 

155.83 

121.79 

155.83 

121.69 


174.5 

165.2 

172.3 

228.5 
228.0 

176.3 
237.9 

188.6 
248.1 
193.7 


Table  A-6.  Data  From  Frank  and  Zook  (1990) 


Test  No. 


(G/cm^) 


L 

(mm) 


1,865 

2,365 

2,409 

2,653 

2,742 

2,746 

3,335 

3,580 

4,415 

4,525 


L/D 

P 

(mm) 

P/L 

V 

(m/s) 

1 

12.5 

0.82 

906 

1 

15.0 

0.98 

1,027 

1 

15.0 

0.98 

1,078 

1 

19.0 

1.25 

1,283 

1 

17.2 

1.13 

1,286 

1 

20.5 

1.34 

1,480 

1 

20.0 

1.31 

1,481 

1 

22.0 

1.44 

1,527 

1 

20.3 

1.33 

1,572 

1 

21.5 

1.41 

1,573 

1 

22.0 

1.44 

1,578 

1 

21.5 

1.41 

1,600 

1 

24.5 

1.61 

1,642 

1 

25.0 

1.64 

1,748 

1 

26.7 

1.75 

1,853 

1 

26.5 

1.74 

1,864 

1 

25.0 

1.64 

1,867 

1 

99.3 

3.90 

4,881 
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